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We investigate tlie pliase diagram of QCD at small chemical potentials, i. e., when chiral and 
flavour symmetry breaking involves the pairing of a quark and an antiquark. The phase diagram 
of two-flavour QCD at small chemical potentials involves chiral symmetry restoration and charged 
pion condensation. We extend previous studies of the topology of the phase diagram, in sections 
with high degree of symmetry, to the physical case of fully broken flavour symmetry, using generic 
thermodynamic arguments. We argue that the extension is unique and present the result. In three 
flavour QCD the phase diagram for chiral symmetry restoration is less well constrained. However, 
we argue that present lattice data allows just two different phase diagrams, which we discuss. 

PACS numbers: 



I. INTRODUCTION 



More than thirty years after the first discussions about a phase transition in QCD |Jj] , only small portions of the phase 
diagram have been explored. Although the complete phase diagram of QCD is of high dimensionality, experiments 
can at best explore a three dimensional section of the full phase diagram. Further, collisions of heavy-ions have only 
a single control parameter, the CM energy, ^/S. As a result they explore a single line in this three-dimensional phase 
1^ . diagram. By varying the nuclei being collided, one could, perhaps, extend the search to a small patch around the 
line. The field is wide open for new ideas on experimental coverage of the QCD phase diagram. 

Theoretical work is no less constrained by the tools of the trade. In regions of high symmetry (for example, when 
the quark masses vanish), universality arguments have been used to put constraints on the phase diagram. Such 
^ ' arguments are realized in models, for example, effective meson models, four-Fermi models or random matrix models, 
\l I which have the same symmetries as QCD. The resulting predictions of universal properties, e. g., the order of the 
£]0 transition and critical indices, are expected to coincide with QCD. Since the locations of phase transitions are not 
■ universal, models should be used to constrain the topology of the phase diagram rather than quantitative predictions 
of the location of phase boundaries or critical points. When the symmetries are broken, as in the real world, the 
usefulness of these models is curtailed further. Weak coupling methods for QCD give precise quantitative predictions, 
but for high temperatures and densities, when the QCD coupling is small enough. Lattice computations were long 
\ confined to the region with vanishing chemical potentials, extensions to finite chemical potential being constrained 
by the fermion-sign problem. The first systematic non-perturbative treatments of QCD at non-vanishing chemical 
potential using lattice methods have now begun, and the first results are now available ^3j|. In spite of these limitations, 
tremendous progress has been made. As we demonstrate in this paper, known results can now be uniquely exended, 
}^ ' using only thermodynamic arguments, to yield the topology of the full phase diagram of two flavour QCD and strongly 
9^ . constrain it for three flavours. 

The topology of the phase diagram of QCD is constrained by its symmetries. It is well-known that QCD possesses 
a set of approximate global symmetries, called flavour symmetries, related to the phases of quark wavefunctions. 
QCD with two flavours of massless quarks would possesses a chiral symmetry SUl{2) x SUr{2) (L and R are 
transformations on left and right handed quarks respectively). The up and down quark masses (m„, nid respectively) 
break this symmetry. Since m = (m„ + md)/2 is non-zero the chiral symmetry is broken to the diagonal vector 
symmetry, S'[/y(2), called isospin. Since m is of the order of a few MeV, and much smaller than the scale, Aqcd, 
chiral symmetry is approximately valid, being broken at the level of 5-10%. The mass difference Am = rud — m„ is 
non-zero but small, thus violating isospin symmetry by a small amount. There is also a semi-light quark flavour, the 
strange, which has mass, nis, comparable to Aqcd- Including this extends the chiral group to 5'J7z,(3) x SUii{3) i|. 
According to data, this symmetry is broken at the 25% level down to the two-flavour symmetry SUl(2) x SUr{2) [5|. 
In addition to these approximate flavour symmetries there is also the exact phase symmetry UsiX) whose charge is 
the baryon number. The axial phase symmetry, Ua{^), is broken at T = by instantons. There is mounting evidence 
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that this symmetry is not restored through a phase transition In this paper we shall assume that there is no 
C/^(l) restoring phase transition in QCD. 

The large global symmetries of QCD can be broken in many ways, thus giving rise to a complicated phase diagram. 
In this paper we confine ourselves to the phase diagram at small chemical potentials. By this we mean that the 
order parameters involve pairing of quarks and antiquarks. At larger chemical potentials there are other interesting 
phases where the pairing could be between two quarks We do not examine these phases in this work. A range of 
intermediate chemical potentials may exist where both quark-antiquark and quark-quark pairings need to be taken 
into account [1]. If this is so, then some details of the phase diagrams presented here would have to be extended. 

Since the phase diagram is structured around the breaking of chiral and flavour symmetry, one might expect that it 
is independent of the number of colours, iVj,. This is correct for > 3. For the specific case of N^, = 2, however, the 
fact that quark representations are real means that the chiral symmetry is enhanced. As a result, the considerations 
below do not apply to Nc = 2. There is a large body of literature on this particular case, and we refer the interested 
reader to a recent review Q. Interestingly enough, as Nc — *■ oo the fact that a baryon contains Nc quarks implies that 
the region of small chemical potentials, in the technical sense adopted here, increases to /i cx Nc- As a result, several 
interesting new phases open up in the hadronic regime and may be studied using different order parameters [lol |. 

We argue here that thermodynamic considerations allow us to extend presently available knowledge to large parts 
of the parameter space of QCD and enable us to build a qualitative picture of the complete phase diagram of QCD for 
small chemical potential. We deduce the topology of the three dimensional slice of the phase diagram of two-flavour 
QCD which may be accessible to experimental tests. We also indicate how these arguments allow us to constrain the 
phase diagram of QCD with up, down and strange quarks. 

The plan of this paper is the following: in the next section we briefly review the Gibbs' phase rule in the form that 
we will use it. The two sections following that deal with Nf — 2 and 3 respectively. The final section contains a 
summary of our results. 



II. THE GIBBS' PHASE RULE 



In this paper we investigate the phase diagram of QCD using an essential tool of thermodynamics: the Gibbs' phase 
rule. Before stating the rule we recall that a system in thermodynamic equilibrium is fully described by a certain 
number of extensive thermodynamic quantities. This set of extensive quantities, among which we must always count 
the entropy S and the energy E, serve as coordinates in the so-called Gibbs space. The thermodynamically stable 
states of a system are in one-to-one correspondence to a convex surface E(S, • • •) in Gibbs space. The thermodynamic 
intensive quantities are derivatives of E with respect to one of the other extensive quantities, the derivative with 
respect to S being T. By the process of taking Legendre transforms of E with respect to each of the intensive 
variables, one reaches a description of thermodynamics in terms of the intensive quantities and a free energy, Q, which 
is extensive. A phase diagram is obtained by projecting out information on Q, and describes the regions in which 
different phases of a system are thermodynamically stable. 

By construction almost all points in the phase diagram correspond to one pure phase. This is the one with the 
lowest free energy among all possible phases. If the system has more than one phase, then at some points in the phase 
diagram two phases may coexist. If we label the phases by a and 6, then the condition for coexistence is the equality 
of free energy densities in the two phases, 

ga{T,fiu,mu, ■■■) = gb{T,fiu,mu, ■ ■ •), (1) 

If the dimension of the phase diagram is D, then two phases coexist along solutions of the above equation, i. e., 
generically along hypersurfaces of 13 — 1 dimensions. Three-phase coexistence requires simultaneous equality of three 
free energies, and hence occurs along hypersurfaces of 13 — 2 dimensions. This surface of 3-phase coexistence is clearly 
the intersection of three surfaces of 2-phase coexistence, obtained by taking the phases pairwise. V phases generically 
coexist along hypersurfaces of £) + 1 — T' dimensions, which are the intersection of P surfaces of two-phase coexistence. 
This is one form of the statement of Gibbs' phase rule. 

A two-phase coexistence surface either has no boundary or ends in a surface of one lower dimension, D — 2, called a 
critical surface. Similarly, a 3-phase coexistence surface may have a boundary. Such a boundary is a D — 3 dimensional 
surface called a tricritical surface. Since the 3-phase coexistence surface is the intersection of three 2-phase coexistence 
surfaces, the tricritical surface is the intersection of three critical surfaces Jjj. Boundaries of 7^-phase coexistence 
surfaces, when they exist, are called P-critical surfaces, and can be viewed as intersections of P-critical surfaces. 
Clearly a P-critical surface is one at which V different phases simultaneously become indistinguishable. Thus another 
form of Gibbs' phase rule states that when the dimension of the phase diagram is D, then there may be 13 — 1 
dimensional surfaces of two-phase coexistence (first order phase transitions), D — 2 dimensional critical surfaces, 
D — 3 dimensional tricritical surfaces, etc.. 
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The most well-known application of the rule is to phase diagrams of chemically pure substances (for example, 
water), characterised by two intensive parameters T and the pressure P. Since D = 2, one has lines of two-phase 
coexistence. Such lines can either end in critical points, or two such lines can meet at isolated triple points (i. e., 
points of three phase coexistence). Both these possibilities are realized in the well-known phase diagram of water. A 
different example is in mixtures of two chemically pure substances (water-alcohol, binary alloys, QCD with quarks, 
etc.), where an additional intensive parameter, the chemical potential /i, makes the phase diagram three dimensional. 
In this case one has surfaces of two-phase coexistence, bounded by critical lines. Three surfaces can intersect along 
lines of three-phase coexistence. The end point of such a line is a tricritical point, and one can view this point also as 
the crossing point of three critical lines. The realization of these possibilities in He'^-He''' mixtures was treated in plj . 

If the phase diagram is known only in some part of the space of intensive variables where calculations are tractable, 
then the Gibbs' phase rule allows us to constrain the possibilities that arise from extrapolations to larger parts of the 
space. As we show in the coming sections, the constraints arc enormous, and sometimes they determine the topology 
of the phase diagram completely. 

In QCD a choice of the extensive quantities can be the order parameters for the breaking of chiral symmetry (one 
condensate for each flavour of quarks and the net number of quarks of each flavour). Conjugate to these are the 
couplings — the quark masses and the chemical potentials respectively. Adding to this set the total energy, one sees 
that for Nf = 2 the phase diagram has 5 dimensions. For Nf — 3 the strange quark mass and chemical potential 
are added, as a result of which the phase diagram is 7 dimensional. Experiments are constrained to work with given 
quark masses, and hence explore a phase diagram of dimension 3 for Nf — 2 and dimension 4 for three flavours. 

This is true if one examines a strongly interacting system on a time scale much larger than the slowest strong 
interaction related relaxation time, but much smaller than the time scale of the flavour changing weak interactions. 
This is the case in heavy-ion collisions, where the expanding fireball cools and freezes out on time scales of the order 
of a few fermis, whereas the typical strangeness changing time scale is many orders of magnitude larger. This might 
seem to indicate that the full 4 dimensional phase diagram can be explored. However, in the initial state the only 
conserved quantum numbers are the baryon number, B, and the electrical charge, Q. Since these two, and the total 
energy, are the only quantities which can be manipulated in experiments, the three dimensional slice of the phase 
diagram conjugate to these variables is the only part of the full phase diagram which is accessible to any possible 
experiment. As is well-known, other parts of the Nf = 3 phase diagram may become accessible in situations at high 
densities and relatively low temperatures, where Fermi-blocking of the light quark states may increase the lifetime 
of strange quarks. This situation may well be realized in compact stellar objects. In the early universe, around the 
time of the QCD phase transition, the baryon density is much higher than now, but still small enough that /is <C T. 
Also, the inverse Hubble time is much smaller than it is now, but much larger than weak the weak interaction time 
scale. Hence, in the early universe, both strong and weak interactions should be considered to be in equilibrium, so 
that the phase diagram reverts to being 3 dimensional. Unfortunately, the chemical potentials in the early universe 
are so small that essentially only the temperature axis is explored. 



III. TWO FLAVOURS 



For Nf — 2 QCD the five intensive quantities can be chosen to be T, fj,dj nriu and m^. A common alternative 
choice is to use /^s, conjugate to the baryon number, B, and /i/, conjugate to the third component of isospin, I^. 
Moving between the ensembles corresponding to these is straightforward; detailed formulae are given in [l2l |. It is 
useful to note that these intensive quantities are conjugate to the extensive thermodynamic variables 

= 1^, s„ . (2) 

and similar relations for down quarks, where Nu is the net number of up quarks and s„ is an extensive quantity whose 
density is the up quark condensate (uu). 



A. Finite temperature 



The best-known part of the phase diagram for Nf = 2 QCD is the temperature axis for m — Am — /.is = l^i ~ 
0. The action has the chiral symmetry SUl{2) x SUfi{2). At T = this is spontaneously broken to the isospin 
part, SUv{2), resulting in three massless pions arising as the Goldstone bosons of the broken symmetry. In the high 
temperature limit the symmetry is restored. The symmetry of the theory is isomorphic to 0(4), which is known to 
have a second order phase transition. Universality then indicates that there is a critical temperature, Tc, at which 
a second order phase transition occurs from the chiral symmetry broken phase to the symmetric phase with critical 
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exponents in the 0(4) symmetry class. The order parameter is the isoscalar chiral condensate, s = (s„ + S£;)/2, which 
changes from a non-zero value at low temperatures to zero at T^. For T < Tc the order parameter is negative for any 
positive m, and flips sign when m in negative. Hence there is a line of first order phase transitions for m = and 
T < Tc, which ends in the critical point at Tc. There is no finite temperature phase transition at finite m. 

Lattice computations are performed at finite quark mass, and are consistent with this picture since they observe 
only a cross over [s^]. Chiral extrapolations have given mixed results: there is no agreement on the critical exponents 
which are seen. This could be a finite lattice spacing artifact. Several studies have used staggered quarks which, for 
m = 0, have only an ?7(1) x U{1) chiral symmetry at finite lattice spacing. The chiral transition in this case should 
lie in the 0{2) universality class. It turns out to be hard to distinguish 0(4) and 0(2) exponents Lattice studies 
are consistent with both [l5|. Studies with Wilson quarks show results consistent with critical exponents in the 0(4) 
universality class [iGj . The situation can be clinched in simulations with dynamical overlap quarks, which realize the 
full chiral symmetry. However, these simulations are currently in their infancy. 
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FIG. 1: The phase diagram in the T-^B-m space [T^] (only the side for m > is shown; a symmetric part of the figure for 
m < is obtained by switching the sign of the order parameter s). A critical line in the 0(4) universality class starts from the 
chiral phase transition, Tc and is the boundary of the first order surface separating the phases with opposite signs of s. The 
tricritical point at (T*,^*) for m — Q lies at the end of the triple line (dashed curve). At the physical quark mass one expects 
a critical end point, (T^,/!^), in the Ising universality class. This lies on the critical line which is the boundary of the first 
order surface separating the phases with non-vanishing and zero value of the order parameter. 



B. The T-HB-m section with and without isospin symmetry 



Recently attention has been focused on the three dimensional section of the phase diagram with varying T, /is 
and m and fixed Am = /i/ = 0. Using the order parameter s one finds the phase diagram shown in Figure [1] The 
plane fj.B = was described in the previous subsection. The first order line and its 0(4) critical end point described 
there extends naturally to three dimensions, becoming a surface of 2-phase coexistence (corresponding to a sign flip 
in s) with a critical line as boundary in the m = plane. For positive m there is a 2-phase coexistence surface across 
which a large negative value of s changes discontinuously to a smaller value. For m < there is a symmetrically 
placed coexistence surface with the sign of s reversed. These three 2-phase coexistence surfaces intersect in a triple 
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line in the m = plane. Since there is a critical line bounding one of these surfaces, it must meet the triple line at a 
tricritical point. We are led to the conclusion that all three 2-phase coexistence surfaces are bounded by critical lines, 
and the tricritical point, at (T*,fi*g), is their common intersection point. 

One fact about this phase diagram is worth emphasizing. The quantity s is an order parameter related to symmetry 
only on the m — plane in Figure [U since it is exactly zero everywhere above the triple-phase and critical lines. At 
finite m, above the coexistence surface it is not exactly zero, but has a value proportional to m. Below this surface 
it has a much larger value. Across the coexistence surface it has a discontinuity, and hence can be used to flag a 
first order phase transition even at finite m. However, along the wing critical lines, s is not an order parameter in 
the sense of being zero in one phase and non-zero in another. This is related to the fact that, for m ^ Q, the phase 
transition is not related to a symmetry. However, along the critical line there are long range correlations between 
local fluctuations in s, or, equivalently, through a fluctuation-dissipation theorem, there are divergent susceptibilities 
of s. This kind of critical behaviour is said to be a liquid-gas transition, and expected to be in the Ising universality 
class. 

The slopes of the coexistence surface can be determined by a generalization of the Clapeyron-Clausius equations. 
For the case at hand we can write 

dGi = SidT + Bid^B + Sidm, (3) 

where the subscript i refers to the two phases (say, a and b) which are at equilibrium along the coexistence surface. 
As one moves along the surface, the changes dT, dfj^s and dm are related by the fact that the free energy density 
is the same in the two phases, eq. ([T]). In order to define a density in a relativistic theory, one notes that all four 
extensive quantities scale similarly, and hence any one of them can be used as a normalization. In determining the 
slope of the coexistence surface, one has to hold one of the three intensive quantities fixed, so one can choose to use 
the extensive quantity conjugate to it to perform the normalization. 

Along lines of constant to, one can then normalize all extensive quantities by s to get the relation 



dfiB 



dT 



J_ Sg/Sg - Sb/Sb 

B l/s,-l/sb 



(4) 



Since the right hand side is independent of the sign of the chiral condensate, one can use the magnitude |sa^6| in the 
expression. Choose a to be the low temperature phase and b to be the high temperature phase. Then, as is shown 
by lattice computations, Sg < Sb and \sg\ > \sb\. As a result, the expression on the right is negative, implying that 
with increasing T the surface of coexistence moves to smaller /i, as shown in Figure [TJ 

Since the ordering of the entropy and chiral condensate determine the slope, the figure is consistent with the idea 
that the low temperature phase contains hadrons, whereas the high temperature phase contains quarks [l^ . At 
finite to one then has a line of first order transitions starting at T = and a large fi b and ending in a critical end 
point at (T^,/i^) which is in the Ising universality class. These arguments, based on the Gibbs' phase rule and the 
Clapeyron-Clausius equation, reproduce the results of earlier model studies ^vf], which spurred recent developments 
in lattice QCD Q and gave rise to experimental interest in the search for the QCD critical point p^ . 

One can extend such arguments to the other two slopes. Along lines of constant ^b one finds 



dT_ 
dm 



(5) 



Given the relative magnitudes of the entropies and chiral condensates, one finds that this slope is negative, i. e., at 
constant fiB, the temperature of coexistence decreases with increasing quark mass. Along lines of constant T one 
finds the slope 



dfiB 



dm 



1 K\/Sg-\Sb\/Sb 

B 1/Sg - 1/Sb ' ^ ' 



which is positive. To the best of our knowledge, arguments about these slopes based on the Clapeyron-Clausius 
equations are new. They imply that with increasing quark mass the critical end point, (T^,/i^), moves to larger T 
and smaller fis, indicated by current lattice data [l9|. 



C. The section T-jiB-f^i when Am — 



The three-dimensional slice through the Nf = 2 phase diagram which is obtained for a fixed non-vanishing value 
of TO and Ato = (including the special slice which also has ^b = 0) has been extensively examined in the literature 
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[20l. [21I [2^. This case has one very attractive feature: the quark determinant is positive definite, and hence lattice 
simulations are possible, and they can be used to check arguments using global symmetries. The drawback is that 
the symmetry which makes the quark determinant positive definite is responsible for making the phase diagram 
non-generic. 

The plane of T-jij, for arbitrary m and vanishing 11 b and Am, was first discussed in ^2QJ . Using an effective theory, it 
was argued that along the line T = 0, there should be a phase transition which results in the development of a charged 
pion condensate, p = ipj^Tiip. It was argued that the phase transition should be critical and in the 0(4) universality 
class with critical fij = ttItt, and that an 0(2) critical line emanate from it. Now, the Gibbs phase rule does not 
generically allow a critical line in a two-dimensional phase diagram. However, subsequent lattice computations [2lj 
saw instead a first order line at finite T connecting to a critical or tricritical point with /Zj oc mT^. 

The phase diagram was extended to the three dimensional slice T-^B-^^I for finite m and Am = by three different 
methods in [23 |. In all these works, critical lines were obtained in the ni-fJ,B surfaces for generic T. Thus there are 
critical surfaces in this three dimensional phase diagram. Such a violation of the Gibbs phase rule may occur in 
subspaces of high symmetry. For example, in the plane m = Am = /i/ = 0, one found the 0(4) critical line shown 
in Figure [H However, these exceptions are non-generic; when extending the phase diagram they do not develop into 
higher dimensional surfaces. This is clear in the above example. In the same way, one is forced to the conclusion that 
the slice Am = is not generic. We consider extensions to Am 7^ next. 



D. Isospin broken by quark masses 



1. The phase diagram in T-fj,i-Am 




FIG. 2: The phase diagram in the T-ni-Am space for generic non-zero m and /ib ~ 0. Since there is a two phase coexistence 
surface on the plane Am — which separates phases with opposite signs of p, this plane also contains a line of triple-phase 
coexistence (dashed, marked Ti). This line may (second and third panels) or may not (first panel) end in a tricritical point at 
T = 0. The last phase diagram is favoured, as we discuss in the text. 



Once isospin is broken by a mass difference between the up and down quarks, an isospin chemical potential no 
longer matches the up quark and down antiquark Fermi surfaces exactly. This breaks the symmetry that kept the 
full quark determinant non-negative, and hence presents the usual problems for lattice simulations. However, it gives 
a generic phase diagram in T-^j. 

Consider the section of the phase diagram for generic non-zero m and fj,B — 0- Changing the sign of Am is 
equivalent to toggling the definition of u and d flavours (in the absence of electroweak interactions) and hence to 
flipping the sign of /i/. Thus, if there is charged pion condensation, i. e., |p| > 0, for some value of /i/, then there is 
a first order transition across Am = 0, in which the sign of p distinguishes the two phases. Thus, along some line in 
the Am = plane one expects triple-phase coexistence, the phases corresponding to p = and two non-zero values 
of p of opposite sign. 

As one goes from small to large /x/ along a line of fixed Am and T, the order parameter p increases. Hence, by 
adapting the argument based on the Clausius-Clapeyron equations for the T-zi^-m phase diagram, one obtains the 
slope drawn in Figure [2l Such a 3- phase coexistence line can be accommodated by the three possible phase diagrams 
shown in Figure O 

The first possibility shown in the figure is generic. The triple-phase line, Ti, starts at T = and continues up. 
There is no critical boundary to the three surfaces of first order transitions. This phase diagram is ruled out by the 
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results of [13, HH which point to some kind of criticahty at T = Am = 0. However, the S'C/y(2) isospin symmetry is 
broken to a vector U(l) when Am ^0. As a result, an 0(4) critical point is not stable under this perturbation. 

The second generic possibility shown in Figure [2] takes this fact into account. The triple-phase line, Ti, ends at a 
tricritical point. Three critical lines emerge from this point. One lies in the symmetry plane Am = and at T 
it is seen to be the 0(4) critical point identified in Along the two wing critical lines the remnant vector U{1) 
symmetry is broken by charged pion condensation. This is closely analogous to the phase diagram in Figure [1] In this 
scenario, as one traverses a line of fixed Am ^ for T = by varying one should go continuously from a phase 
with small p to one with a large value of this condensate. However, at T = and small fij it is not energetically 
feasible for pions to condense, hence there must be a range of fij over which p vanishes. In that case a phase transition 
must occur at T = for any ^j. 

This brings us to the final possibility: the triple-phase line ends at T = in a tricritical point. Two wing critical 
lines emerge from this in the T = plane. The third critical line, which is the 0(4) line of the second scenario, is 
squashed down to a point. In this scenario there is always a phase transition between the charged pion condensed and 
uncondensed phases. This is the favoured phase diagram, since the lattice computation of [2l|, which sees a critical 
or tricritical point at T = and Am = 0, prefers this over the first scenario. It would be good to check this through 
another computation at smaller lattice spacing. 

2. The T-^B-Am phase diagram 




FIG. 3: Two possible topologies for the phase diagram in T-HB-Am for fii — 0. The second one is preferred (see the text for 
details) . 

One might expect that when the isospin symmetry is broken, the flavour singlet order parameter should be replaced 
by the pair s„ and [2^ . If so, the T-jib phase diagram would be extended to the first topology shown in Figure [3l 
This contains a line of 4-phase coexistence, labeled Q in the diagram with an end point which is a tetracritical point. 
In a three dimensional phase diagram one does not generically expect either a 4-phase coexistence line or a tetracritical 
point. We have argued before that non-generic situations might be obtained on sections of high symmetry. However, 
having two different non-vanishing masses breaks the flavour symmetry maximally, so it is impossible to find a more 
generic situation. 

This no-go argument can be evaded by noting that and Am break the fiavour symmetry to the same subgroup. 
As a result, it is possible that the topology is the same in the two sections T-/iB-Am for jij ~ and T-fiB-fJ-i for 
Am = separately. In other words, the 4-phase coexistence line and the tetracritical point lie in the intersection of 
these two sections in this extended four dimensional phase diagram r-//B-/i/-Am. By the same token, such a topology 
is ruled out as non-generic in the T-/zb-^/ section for Am = unless it is also seen in this section. 

Arguing differently, one might expect that for Am/m < 1 the effect of isospin symmetry breaking through the 
quark masses should not change the phase diagram qualitatively. In that case the two order parameters s^, and 



8 



would be redundant, and the two sheets of first order transitions merge into a single sheet of first order transitions 
found with s. This possibility is shown as the second topology in Figure [H 

If this second argument is correct, then one must face the question of what happens at large Am/m. Exactly 
this question was answered at T = using an effective theory in \23\. At finite Am, the explicit breaking of isospin 
symmetry to vector U(l) allows the tt*^ to mix with isoscalars (this allows a change in its mass without changing the 
masses of tt^). When Am/m is large, a neutral pion condensate can form, and, at T = 0, there is a second order 
transition between phases with and without this condensate. In the real world a mass difference between the charged 
and neutral pions is observed, so a mixing of this sort could occur. However, a neutral-pion condensate is not observed 
at r = 0, so one may conclude that the physical value of Am/m is not large enough for isospin breaking effects to 
qualitatively change the phase diagram shown in the second part of Figure [31 

The mixing of neutral pions with isoscalars is a physical effect which does not exist in the models used in 22]. This 
leads us to prefer the second phase diagram shown in Figure [31 The only lattice study to date of thermodynamics 
at finite Am [2^ saw no evidence of two separate crossovers. A computation with a different lattice formulation of 
quarks will be an interesting and useful additional piece of evidence. 

E. The physical phase diagram 



T 




FIG. 4: The phase diagram in the T-jiB-yii space for generic non-zero m and Am is organized by the two order parameters s 
and p. It contains a line of 3-phase coexistence (labeled TS) and two tricritical points (labeled A and A'). 

The physical phase diagram of QCD is the one in which all possible experimentally tunable parameters are shown 
for fixed physical values of the remaining parameters. This is the phase diagram in the space of T-^b-^J-i for generic 
small m and Am < m. After the analysis of the previous subsections, deducing the topology of this phase diagram 
is straightforward. The various phases are distinguished by the order parameters s and p. We discuss the topology 
of the phase diagram (shown in Figure [4]) only in the quadrant with /is./ > 0. The remainder can be constructed by 
symmetry. 
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In the /z/ = plane there is a first order line ending in the QCD critical point. Across this line s changes 
discontinuously. Following the discussion of the previous subsection, we conclude that this line develops into a surface 
of 2-phase coexistence. The QCD critical end point develops into a critical line (in the liquid-gas universality class) 
which is the boundary of this surface. In the /is = plane there is a first order line separating phases with vanishing 
and finite values of p with a liquid-gas critical point at T = 0. This develops into a 2-phase coexistence surface. 

The two first order surfaces must meet along some line, TS. Hence T3 is a line of 3-phase coexistence, and another 
surface of 2-phase coexistence must emanate from it, separating the phase in which both the condensates vanish from 
that in which neither does. The two ends of T3, labeled A and A' , must both be tricritical points. At A the first 
order surface of discontinuity in p must have a boundary. Since the system has no particular symmetry at A, all 
three critical lines meeting at this point must be in the universality class of a liquid-gas transition. At A' there seem 
to be two critical lines, both in the liquid-gas universality class. The third critical line is degenerate, as in the third 
topology shown in Figure [21 

The derivation of the Clapeyron-Clausius equations for the physical phase diagram needs little further comment. 
Along planes of constant /is one can use eq. (jlj, since m and Am are also constant along this line. Along the 
planes where /i/ are constant, the modification of eq. ^ was already discussed in Section fill D H and leads to the 
opposite slope of the coexistence surface for pion condensation. The remaining part involves scaling the expression 
dGi ~ SidT + BdfiB + IsdfJ-i by Si in the phase i, and equating the two free energy densities so defined at each 
coexistence surface. The resulting Clapeyron-Clausius equation is 



The first order surface of chiral symmetry restoration, which is in the phase without a pion condensate, then has only 
weak dependence on /i/. Similarly, the first order surface of pion condensation has weak dependence on /is because 
the baryon density is low in the phase with chiral symmetry breaking. The surface separating the phase with both 
condensates vanishing from the phase where both condensates are present has more complicated behaviour. One 
expects that in the phase without condensates the number densities scale with chemical potentials as -B cx fi^ and 
/3 cx /ij when T is small. In that case this first order surface is roughly linear in and /i/. Lattice computations 
at Am = can substantially improve our knowledge of the functional form of B{T, fiB, l^i) and I^{T, ^b, ^^I)^ and 
therefore the shape of the physical phase diagram through continuation to Am ^ 0. 

The position oi A and A' are constrained by these considerations. It is estimated that the QCD critical end point at 
Am = /X/ = is at fig/T^ ~ 1 with ~ 170-190 MeV. Since the effect of isospin breaking is small, one expects that 
at realistic values of Am the critical end point does not shift drastically. One also knows that at /is = Am = 0, the 
critical point in fij ~ m.^. The corrections at realistic Am are not large. One expects that ~ ~ ~ 100-200 
MeV. Also, from this topology one could deduce that in heavy-ion collisions where /x/ <C m7r one should expect to see 
the QCD critical line close to the estimates of and . 



For Nf — 3 QCD, the intensive quantities at the phase transition are the temperature T, the three quark masses, 
ruu, md and mg, and the three chemical potentials /x«, /id and /is- As before, we can take linear combinations 
of these variables to parametrize the seven-dimensional phase diagram. We will use the mean light quark mass 
m = (m„ -|- md)/2 and the mass splitting between these two. Am = m„ — rud- It is useful to construct the baryon 
chemical potential, ^b and two other combinations, as outlined in [l^ . 

The three fiavour phase diagram has barely been explored, and we shall not be able to do justice to the many kinds 
of pairings that it may possess. In this first work we restrict attention to a question which has been asked recently: 
is the two-flavour phase diagram for chiral symmetry restoration (Figure [1]) a good guide to the corresponding phase 
diagram in QCD with additional strange quarks? 



Universality arguments [3] lead us to expect that the chiral phase transition in three flavour QCD at zero chemical 
potentials is of first order. The order parameter used is the three-fiavour condensate. In the three-fiavour chiral 
limit, since the action is blind to flavour, any linear combination of three flavoured condensates is as good an order 
parameter as any other, so one could as well continue to use the two-flavour order parameter s to study the phase 
structure. When the strange quark becomes massive with the light quarks remaining massless, one does not expect 




B 



(7) 



IV. THREE FLAVOUR 



A. Vanishing chemical potentials 
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FIG. 5: The phase diagram of QCD for vanishing chemical potentials and Am = (on the left) as observed with the light 
quark condensate. In different parts of this phase diagram different flavour symmetries are obtained: the diagonal plane has 
exact 3-fiavour symmetry, the "back" plane has Nf = 2, the far corner is the quenched theory, and the plane at the far right 
has Nf — 1. On the right is the corresponding flag diagram, obtained by projecting the phase diagram to the plane of T = 0. 



symmetry restoration in the strange sector while the symmetry remains broken in the light quark sector. Similar 
expectations hold whenever the strange quark mass is heavier than the up and down quark masses. Hence, in this 
whole region one expects that the same phase structure is seen whether one uses the three flavour condensate or s. 

The phase diagram of QCD with three quarks has been explored in several lattice computations [1^, [H, [23, HI] 
where the two light quarks are degenerate, i. e.. Am = 0. In this case the phase diagram in the section m-ms-T is of 
the kind presented in Figure [3 In these studies there is only one chiral phase transition observed, and not separate 
ones for strange and light flavours, in conformity with the above discussion. For infinite and m — 0, one has two 
flavour chiral symmetry, and hence a second order transition in the 0(4) universality class. When nis = m, i. e., three 
flavour symmetry is exact, one has a line of flrst order transitions up to a critical end point, which lies in the Ising 
universality class [23|. Moving off the exact Nf — 3 line (termed going to Nf = 2 + 1 in the literature), this develops 
into a surface of first order transitions bounded by a critical line. Since this line is in the liquid-gas universality class 
[131, it must meet the 0(4) line in a tricritical point. The nature of this critical line is crucial in establishing a phase 
diagram. If it is in the liquid-gas class, as expected, then this would be visible with any formulation of lattice quarks. 

The phase diagram shown in Figure [5] can be probed entirely with the light quark condensate. There are, of course, 
other phase transitions in this phase diagram — for example the deconfining transition in the corner which contains 
the quenched theory. This is observed with a different order parameter, the Polyakov loop, and hence does not appear 
in the phase diagram shown. 

These results are often presented in a figure (the second of Figure [5]) which is the phase diagram projected down 
to the plane m-ms- This is not a phase diagram since a point on it does not represent a stable thermodynamic 
phase. Instead, each point on this diagram represents whether or not there is a phase transition at some T "above" 
it. Unfortunately there seems to be no name for such an useful diagram, and we are forced to invent a name for it. 
In the rest of this paper we call it the "fiag diagram" , because each point flags whether there is a phase transition 
above it, and if so, the order of the transition. 

One of the main open questions about the flag diagram is the location of the physical point. All computations 
indicate that this lies deep in the crossover region [2^ [2^ . A second important question is about the location of the 
tricritical point. What is the value of m^'^? If the physical strange quark mass is less than mf'^ then the Nf = 2 
phase diagram will be a good guide to the phase diagram of real QCD. If, on the other hand, the physical strange 
quark mass is larger, then the phase diagram in the real world could be significantly different. One naively expects 
that ml'^ ~ Aqc£,; since the physical value of TOs is also similar, this is an interesting and wholly non-perturbative 
question. 

One estimate of m^'^ that we are aware of was made using a linear sigma model with parameters fixed by hadronic 
data iSOj. This finds that the K meson mass at the tricritical point is at least 1700 MeV, far in excess of the physical 
value. While this would seem to indicate that the physical strange quark mass is lighter than m^'^, there are indications 
that further work may be needed in order to pin down m^'^. For example, in [30] it is estimated that the critical end 
point along the Nf = 3 fine occurs at — = 110 ± 20 MeV. The lattice computations of [l^l indicate that the 
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FIG. 6: Possible flag diagrams for three flavour QCD are organized around the shape of the tricritical line (dash-dotted line). 
The surface of 0(4) criticality which lies in the plane with m = meets a surface of Ising criticality along this line. The two 
possibilities shown are distinguished by the approach of the tricritical line to the hb = plane. The second possibility includes 
the case where the tricritical line goes to infinite somewhere between its two fixed ends. 

critical end point may be at = — 192 ± 25 MeV. Whether the discrepancy is due to an insufficiency in the 
model of ISOll or the multiple extrapolations performed in [27| . or a combination of the two, remains to be tested in 
future. In [3l| the estimate m^'^ ~ 500 MeV is presented. Indications that the cutoff effects at comparable cutoffs are 
large come from the results of [2^, [2^ . 



T 




FIG. 7: The phase diagram for fixed m*"""" < < ml'^ and Am = /ii = jiy = in the second scenario shown in Figure |5] 
The phase diagram for fixed m is shown in the inset. 
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B. Finite /is 

Another slice that is interesting from the point of view of forthcoming planned experiments is the four-dimensional 
T-^B-m-ms slice obtained with Am — fii — /ly = 0. This was investigated recently on the lattice [2^.[3]|. and a flag 
diagram obtained by projecting out the T axis was used in the presentation. These two studies differ in the conclusions 
reached: one [l^ favours the first flag diagram shown in Figure [6] the other [3l| , we argue below, is consistent with 
the second possibility shown in the same figure. These are computations which present many technical challenges. 
We remark later on one way to resolve the issue. 

Two sections of this phase diagram have been discussed — the section with /.ib = was discussed in the previous 
subsection, and the section with nis — was discussed in the section on the Nj — 2 phase diagram in T-^B-m- Both 
contain a line of 0(4) critical points, and a line of Ising critical points, meeting at a tricritical point. The two lines 
of 0(4) critical points are in fact the same, since they both originate in the finite-T transition for the N f = 2 chiral 
symmetry restoration. In the 4-dimensional phase diagram one should, therefore, see a single critical surface of 0(4) 
transitions. The boundary of this critical surface is the tricritical line, the two ends of which show up as the two 
tricritical points we have already discussed. At this tricritical line the 0(4) critical surface is glued on to surfaces of 
Ising criticality. The tricritical line can also be viewed as the boundary of a surface of triple phase coexistence. This 
much is a direct consequence of the Gibbs phase rule. 

The Gibbs phase rule also allows a tetracritical point in a generic four dimensional phase diagram. A tetracritical 
point is the end point of a line of four-phase coexistence. With a single scalar order parameter, such as s, it is not 
possible to distinguish four phases. Therefore, we conclude that a tetracritical point should not be seen in the four 
dimensional phase diagram observed with s. 

There remain two generic possibilities for the shape of the phase diagram and they are shown in terms of the 
corresponding flag diagram in Figure [SI The first possibility is that the tricritical line lies entirely on the side 
> ra\'^ . The second is that the tricritical line continues to some m'"''" < to^"^ and then bends back towards 
m^'^. In the first scenario, the critical line moves to larger m and as /is is increased by a small amount from 
/xb = 0. In the second, the opposite movement occurs until a large value of /is. This is precisely the observation 
which distinguishes the two lattice studies in [2^ [3l| . It would be interesting to push lattice studies to smaller cutoffs 
and decide which scenario is actually obtained in QCD. 

In the first case the phase diagram for > m^'^ is qualitatively like the N f = 2 phase diagram, and for rus < rn^J^, 
similar to Nf = 3. In the second case for m^"^ > TOj > to*'"'" there are two sheets to the phase diagram at fixed nis' 
one similar to the Nf — 2 phase diagram, the other like the N f = 3 phase diagram, with the two tricritical points 
joined by an 0(4) critical line, as shown in Figure [T] The inset in the figure shows the unusual phase diagram for 
fixed non-zero value of m. It would be very interesting if QCD is eventually seen to yield such a phase diagram. 

V. CONCLUSIONS 

The phase diagram of QCD is complicated. The large global symmetries of QCD imply that there are many 
thermodynamical intensive parameters, i. e., free parameters which enter the Lagrangian. As a result, there are many 
different kinds of pairings of quarks that can arise as these parameters are changed. In this paper we have considered 
only small chemical potentials, defined to be the region where the relevant pairings are between a quark and an 
antiquark. 

The topology of the phase diagram for QCD with two flavours of quarks has been investigated earlier by many 
authors, mainly in sections of partial symmetry. In this paper we have extended these considerations to the case of 
physical quark masses, which break both the chiral symmetry and isospin symmetry. Explicit breaking of isospin 
symmetry by quark masses leads to mixing of isovectors with isoscalars and hence influences the phase diagram 
(see Section fill D 2}i . This piece of physics is incorporated into the phase diagram that we construct. The topology 
of this physical phase diagram is strongly constrained by thermodynamics — essentially the Gibbs phase rule (see 
Section [n|. Our results are given in Section fill El and summarized in the phase diagram shown in (see Figure [4]). 
The physical phase diagram contains a triple line at which phases with and without the chiral and pion condensates 
coexist. Deducing actual values of the critical and tricritical points in this diagram is a much harder task. However, 
appropriate Clapeyron-Clausius equations along with the results of current lattice computations place various 
constraints on them. Improving these constraints through future lattice computations is possible without extreme 
effort. 

The phase diagram of three flavour QCD could be more complicated, since there are even larger global symmetries. 
Strange quarks can decay into the light quarks through the weak interactions, so one may ask whether this phase 
diagram is relevant to any physical system. In the context of heavy-ion collisions, where the system is expanding 
so fast that weak interactions do not come into play in the thermodynamics, strange quarks are indeed relevant. 
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Further, in the physics of compact stars (involving large chemical potentials), where weak decays of strange quarks 
are suppressed due to Fermi-blocking of the light quark states into which they could have decayed, strange quarks 
also become relevant to thermodynamics. Thus the study of the three flavour phase diagram is interesting. 

Comparatively little is known at present about this phase diagram. Here we restrict attention to the breaking and 
restoration of chiral symmetry. The Gibbs phase rule places strong constraints on the topology of the phase diagram. 
As the strange quark mass is increased keeping the light quarks massless, there is a tricritical value, rTi^"^, at which the 
finite temperature chiral phase transition changes from being first order to second order (for rus > m^'^). Presently 
available observations on the lattice [29, 31] can be patched together into two possible phase diagrams which we have 
discussed in Section IIV Bl The three-flavour phase diagram can be improved substantially with lattice computations 
in the near future. 

Apart from the results on the QCD phase diagram we believe that this paper demonstrates an useful technique. 
Usually phase diagrams have been examined by writing down a Ginzburg-Landau theory. As remarked before, these 
give useful predictions for universal quantities. However, when all symmetries are broken, the usefulness of such a 
theory is severely curtailed. In such situations, we found that the Gibbs' phase rule is a handy tool for exploring the 
topology of phase diagrams. When supplemented with the Clapeyron-Clausius equations, it may even be possible to 
make some quantitative statements of the kind that have been explored in this work. 

It is a pleasure to thank Mike Creutz, Saumen Datta, Philippe de Forcrand, Rajiv Gavai, Bengt Petersson, and 
Misha Stephanov for discussions. This work was supported by the Indian Lattice Gauge Theory Initiative. 
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